arXiv:1506.05699v2 [physics.soc-ph] 11 Sep 2015 


Stochastic Evolving Differential Games Toward a 
Systems Theory of Behavioral Social Dynamics 

Giulia Ajmone Marsan 

Organization for Economic Co-Operation and Development 
Paris, France 

giulia.AJMONEMARSANSoecd.org 

Nicola Bellomo 

Department of Mathematics, King Abdulaziz University 
Jeddah, Saudi Arabia 

nicola.bellomoOpolito.it 

Livio Gibelli 

Department of Mathematical Sciences, Politecnico di Torino 
Torino, Italy 

livio.gibelliSpolito.it 


September 14, 2015 


Abstract 

This paper proposes a systems approach to social sciences based on 
mathematical framework derived from a generalization of the mathemati¬ 
cal kinetic theory and on theoretical tools of game theory. Social systems 
are modeled as a living evolutionary ensemble composed by many in¬ 
dividuals, who express specific strategies, cooperate, compete and might 
aggregate into groups which pursue a common interest. A critical analysis 
on the complexity features of social system is developed and a differential 
structure is derived to provide a general framework toward modeling. 
Keywords: Kinetic theory, active particles, stochastic differential games, 
evolution, learning. 

AMS Subject Classification: 82D99, 91D10. 


1 Motivations and Plan of the Paper 

This essay aims at developing a modeling approach of social agents viewed 
as behavioral, evolutionary, complex systems. It provides a critical overview 
on the contribution of mathematical sciences to the modeling, qualitative and 
computational analysis of social systems with some implications also in the field 
of economics. Subsequently, a systems approach is developed and applied to 
selected case studies. 
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Various hints coming from a radical change observed in social and economic 
sciences [16] suggest to take into account the complexity features belonging to 
socio-economic systems MUITllIlillliiHlIMiniT]. The need of including 
in the modeling approach the specific features of living system is well motivated 
in the book by Baley [18] . where new concepts of entropy and equilibrium are 
developed. The same author, in a second book |19j proposes several new ideas 
toward a theoretical synthesis for a systems approach to social systems. An 
additional field of application includes a systems approach to urban planning, 
as an example papers [231 [231IS3] , propose a systems approach of cities viewed 
as complex systems. 

Motivations are mainly induced by a radical philosophical change that is un¬ 
folding in social and economic disciplines. Roughly speaking, the new emerging 
point of view is characterized by an interplay among Economics, Psychology, 
and Sociology, which is no longer grounded on the traditional assumption of ra¬ 
tional socio-economic behaviors. The rationale for that approach, namely that 
Economics can be highly affected by individual (rational or irrational) behaviors, 
reactions, and interactions is widely accepted. 

Therefore, an important hallmark in the modeling of socio-economic systems 
consists in understanding the complexity features of living systems and of so¬ 
cial interactions in particular. Subsequently, mathematical methods need to be 
developed which are suitable, as far as it is possible, to capture such features. 

The recent literature shows that behavioral features can be taken into ac¬ 
count by suitable developments of methods of statistical mechanics [91] , kinetic 
theory evolutionary game theory and applications [26l [271E?] EH] |94l 11121 
ITT^fTMrm] . mean field games |105j . and more specifically by the so-called 
kinetic theory of active particles [37], called KTAP, which combines, as we shall 
see, methods of the classical kinetic theory and theoretical tools of game theory. 

Speculations and research activity in the field of physics cannot be neglected 
by mathematicians, hence the book by Galam m, as well as various papers of 
the same authors such as [7E1IE2], are an important reference for the approach 
proposed in this paper. 

Let us now focus with more detail to evolutionary game dynamics which pro¬ 
vides an important conceptual contribution to this present paper. This topic 
was presented and critically analyzed in the survey paper [94] . where the au¬ 
thors put in evidence the substantial difference compared to the classical game 
theory [nnidn]. Namely rather than players involved in a game who attempt 
to maximize their payoff, a whole population is entangled to pursue an indi¬ 
vidual or collective wellbeing. This concept is well expressed by the following 
quotation extracted from [94] : 

Evolutionary game theory deals with entire population of players, all 
programmed to use the same strategy (or type of behavior). Strate¬ 
gies with higher payoff will spread within the population (this can 
be achieved by learning, by copying or inheriting strategies, or even 
by infection). The payoffs depend on the actions of the co-players 
and hence on the frequencies of the strategies within the population. 
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Since these frequencies change according to the payoffs, this leads to 
a feedback loop. The dynamics of this feedback loop is the object of 
evolutionary game theory. 

Moreover, the authors also observe in that one of the great contributions 
of Mayr is the replacement of typological thinking by population thinking |107j . 
This concept linked to the dynamics of mutations and selection can explain 
various aspect of the theory of evolution [108] . The authors of this present 
paper think that tools of game theory need including, as we shall see, features 
of evolution dynamics, learning processes. An additional important topic, which 
will be widely treated in the next sections, is the role of nonlinear interactions. 

The aforementioned kinetic theory for active particles combines in a gen¬ 
eral mathematical framework all above concepts including the complex role of 
heterogeneous behaviors over the collective dynamics which can include muta¬ 
tions and selection. Our paper focuses on this approaches and the presentation 
closely follows the style of a recent book chapter nni, where five key questions 
are posed by the authors to themselves, as well as to the attention of the reader. 
The answer to such questions contributes to achieve the aim of developing a 
mathematical approach to behavioral social dynamics. 

The presentation of this paper accounts for the knowledge that the research 
activity in the field is still in a preliminary stage and looks for a commonly 
shared strategy. Hopefully, different efforts, including those proposed in this 
paper, might converge to the common objective of developing a systems theory 
approach. 

Bearing all above in mind, the following five key questions are proposed: 

1. How mathematicians should look at a behavioral social dynamics 
and to related new research hallmarks? 

2. Should the derivation of models refer to a general mathematical 
structure suitable to offer a conceptual framework? 

3. How far the existing literature has exploited a “mathematical 
structure” ? 

4. Can a general systems theory be developed and how can it con¬ 
tribute to the interpretation of the so-called “big-data” ? 

5. How and how far, mathematical models can be validated? 

The answer to these key questions is given in the next sections, where math¬ 
ematical tools refer to the KTAP’s methods reviewed and revisited in m- More 
in detail, the contents of this paper is guided by the authors’ opinion that suit¬ 
able mathematical methods can significantly contribute to a deeper understand¬ 
ing of the relationships between individual behaviors and the collective social 
outcomes they spontaneously generate. Such mathematical tools should expand 
the concept of bounded rationality introduced by Simon |12511126] , according 
to the idea that choices made by individuals are not perfectly rational as their 
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rationality is limited by the information they have, the cognitive limitations of 
their minds and the time available to make the decision. Decision-makers in 
this view act as satisfiers who can only seek a satisfactory solution, lacking the 
ability and resources to arrive at the optimal one, see also [53 1130] . 

The alternative approach, which this paper refers to, looks at large social 
systems whose dynamics is based on individual behaviors, by which single indi¬ 
viduals express, either consciously or unconsciously, their own strategy, which 
is often based not only on their individual purposes but also on those they 
attribute to other agents. 

The contents of each section are not limited to technical answers, and a 
broader analysis is given looking ahead to perspectives. Adding further indica¬ 
tions on the contents: 

Section 2 provides, by answering to the first key question, a phenomenolog¬ 
ical interpretation of social systems which exhibit a behavioral dynamics. The 
hallmark, which is followed, is that the modeling approach should be referred 
to a mathematical structure suitable to retain the complexity features of social 
systems. 

Section 3 focuses on the second key question and derives a general structure 
which retains the complexity features studied in the preceding section. Such a 
structure is deemed to provide the framework for the derivation of models. 

Section 4 faces the third key question and proposes a critical analysis of the 
existing literature, also with the aim to show that some applications fall into 
the general structure derived in the preceding section. 

Section 5 proposes a systems sociology approach, where the mathematical 
structure can be adapted to the study of a variety of case studies corresponding 
to different fields of social sciences. Moreover, computational methods and 
the interpretation of large database are critically analyzed. Hence, this section 
specifically refers to the fourth key question. 

Section 6 presents the analysis of a case study originally developed in [32] 
which models the interaction among citizens, police forces, and criminality. It 
is shown how various concepts of the preceding sections can be implemented at 
a practical level. 

Section 7 tackles, focusing on the fifth key question, the delicate problem of 
model validation based on the key problems proposed in the book [^, which 
will be reported in the next section. Validation cannot be fully achieved for 
living systems as these are subject to large deviations. However some useful 
hallmarks can be defined and properly applied. 

Section 8 outlines some research perspectives selected according to the au¬ 
thors’ bias as one of the aims of this paper consists in addressing future research 
activity in the challenging field of mathematics interacting with social sciences. 

2 Toward a Behavioral Social Dynamics 

Although, a model is only an approximation of reality, the contributions of 
mathematics first to physics and subsequently to applied sciences is well rec- 
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ognized and is in continuous evolution to achieve a deeper understanding of 
systems of the real world. 

Mathematics has developed rather sophisticated qualitative and quantitative 
tools for studying the inert matter, where causality principles can be generally 
applied. On the other hand, the modeling of living matter cannot rely on direct 
cause-effect links, because the active ability of individuals to develop behavioral 
strategies and to adapt them to the context, makes observable effects come out 
from often non-evident causes. 

The most concrete example is in the field of biology, where despite a long 
lasting interaction and some successful result, mathematics is still very far from 
understanding the complexity of biological phenomena. The lack of invariance 
principles has been well highlighted in biology [93l 110611108] , while similar ar¬ 
guments can be proposed focusing on the interplay between Mathematics and 
Socio-Economic Sciences. 

This issue is considered in the first key question concerning the conceptual 
interactions that a mathematician should develop in the study of social systems. 

How mathematicians should look at a behavioral social dynamics 
and to related new research hallmarks? 

Mathematics needs transferring what is observed through a phenomenologi¬ 
cal interpretation of social reality into equations, generally of differential nature. 
Therefore it is important selecting a number of key features toward the aforesaid 
mathematical purpose. 

In general, living systems exhibit various complexity features. In particular, 
interactions among individuals need not have an additive linear character. As 
a consequence, the global impact of a given number of entities {field entities) 
over a single one {test entity) cannot be assumed to consist merely in the linear 
superposition of the actions exerted individually by single field entities. This 
nonlinear feature represents a conceptual difficulty to the derivation of mathe¬ 
matical models. More generally, the new point of view presents social science as 
a discipline, which should deal with evolving complex system, where interactions 
among heterogeneous individuals and the interplay among different dynamics 
can even produce unpredictable emerging outcomes [MISS]- 

After these preliminary statements, the following complexity features are 
selected: 

• Ability to express a strategy: Living entities are able to develop specific 
strategies to fulfill their goals depending on their own state and on that 
of the entities in the surrounding environment. Living systems typically 
operate to achieve their well-being. The said strategies are not defined 
once for all, but are modified according to the changing environment where 
individuals operate. 

• Heterogeneity: The ability to express a strategy is heterogeneously dis¬ 
tributed. Ln the case of social (and economic) systems, heterogeneous be¬ 
haviors can play an important role in determining the overall collective 
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dynamics whenever irrational behaviors of a few entities can induce 
large deviations from the usual dynamics observed in rationality-driven 
situations. 

• Nonlinear interactions: Interactions are nonlinearly additive and in¬ 
volve immediate neighbors, but in some cases also distant individuals due 
to the ability living systems to communicate. In some cases, the topologi¬ 
cal distribution of a fixed number of neighbors can play a prominent role 
in the development of strategies and interactions. Self organization is an 
important consequence of interactions. 

• Learning and adaptation: Living systems have the ability of learning 
from past experience. Therefore, their strategic ability and the character¬ 
istics of mutual interactions evolve in time due to inputs received from 
outside. As a result, continuous adaptation to the changing-in-time envi¬ 
ronmental conditions occurs. 

• Selection and evolution: Aggregation of individuals can be viewed, as 
the formation of groups sharing the same interest which may generate 
new groups more suited to an evolving social and economical environment 
as well as groups less suited to the social context. In the course of time, 
some of such groups may disappear, and new ones may become dominant. 

According to these features, it is apparent that the dynamics of a few entities 
does not lead straightforwardly to that of the whole system, because the latter 
manifests, at a collective level, the nonlinear effects of individual interactions. 
This observation can be further understood based on the following statements, 
extracted from [^: 

Living entities express, rationally or irrationally, a certain strategy 
for achieving their own wellbeing. These systems are often complex: 

At large scale collective behaviors appear, which are apparently coor¬ 
dinated but are actually self-organized. Namely, they emerge spon¬ 
taneously without the action of any external organizing principle. 

Sudden deviations from the usual standards may occur, leading to 
even highly unpredictable events with dramatic collective consequences, 
as well demonstrated nowadays by recent facts in our societies. One 
of such events, which is of paramount interest in Social Sciences 
in general, and in Economics in particular, is the so-called Black 
Swan, namely an extreme event, largely unpredictable at a collective 
level, originating from apparently rational and controlled individual 
behaviors fl2ff\/ . 

Remark 2.1 Complexity in living systems is induced by a large number of vari¬ 
ables, which are needed to describe their overall state. Therefore, the number of 
equations needed in the modeling approach might be too large to be practically 
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treated. Appropriate methods need thus to be developed with the aim of reducing 
this structural complexity. At the same time these methods should be multiscale 
in nature since using only a single observation and representation scale is not 
sufficient to describe the dynamics of living systems. For instance, the micro¬ 
scale of individuals affects the scale corresponding to a group of interest, which 
in turn has a direct influence on the dynamics of large systems made up of the 
union several interacting groups of interest. 

Remark 2.2 Living systems composed by many individuals show collective emerg¬ 
ing behaviors that are not directly related to the dynamics of a few entities, but 
are often generated by a kind of swarming intelligence that involves all the in¬ 
teracting individuals 11381 \5‘A At a qualitative level, emerging behaviors are 
often reproduced under suitable input conditions, though quantitative matches 
with the observations are rarely obtained. In fact, small changes in the input 
conditions often generate large deviations. In some cases, these break out the 
macroscopic (qualitative) characteristics of the dynamics, whence substantial 
modifications can be observed. Heterogeneity of individual strategies, learning 
ability, and interactions with the outer environment largely influence such phe¬ 
nomena. According to /5'6] /. these deviations can be interpreted as the requisites 
for a “Black Swan” to develop \129f . More precisely, the following definition 
can be extracted from the cited book: 

A Black Swan is a highly improbable event with three principal 
characteristics: it is unpredictable; it carries a massive impact; and, 
after the fact, we concoct an explanation that makes it appear less 
random, and more predictable, than it was. 

Remark 2.3 Mathematical models can hopefully offer a powerful tool to show 
the evolution of different possible scenarios in a socio-economic systems. They 
might even lead to better understand what possible directions emerging collec¬ 
tive actions may take, depending on different initial conditions and individual 
behaviors. 

All the remarks above have been presented in agreement with the authors’ 
bias. However, it is not claimed that this section offers an exhaustive framework. 
Therefore, we wish to present the key problems proposed by Bonacich and 
Lu [53], with the aim of getting back to them for further speculations after 
having presented the mathematical tools proposed by our paper. 

1. Regarding networks and homophyly birds of feather flock together, but peo¬ 
ple are influenced by those they like. Both these processes result in the 
same outcome, but there is no standard accepted way of separating these 
two processes. 

2. There are lots of models that show how groups arrive at consensus but no 
generally accepted model of how groups become more and more different 
and possibly hostile. Moreover, what are the most important mechanism 
for bringing about cooperations in groups? 
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3. We know that people are affected by their positions in networks, but we do 
not have a variety of models of how people create their networks. We also 
do not have good models for network change and evolution. 

4. There are a variety of measures of centrality in networks, but there are no 
well-established criteria for when one measure is preferable to another. 

5. Are human groups unpredictable because humans themselves are complex 
organisms or because there are truly chaotic dynamics in groups? 

We have reduced the “key problems” from the six proposed in [53] to five, 
simply by grouping the previous second and fourth problem. Actually, the 
authors refer in [^ to the celebrated Hilbert’s problems, which have highly 
engaged the mathematical research activity in the past years and are still object 
of challenging speculations. Although the reference to Hilbert appears to us 
somehow excessive, still we think that these problems offer possible way to 
validate models and theories or, at least, critically analyze them. Therefore, the 
whole matter will be treated again in Section 5. 

3 On the Search of a Mathematical Structure 

This section provides an answer to the second key question, which focuses on 
the search of structure suitable to capture the specific features of social systems 
proposed in Section 2. 

Should the derivation of models refer to a general mathematical 
structure suitable to offer a conceptual framework? 

The answer is given through four subsections. Subsection 3.1 provides a gen¬ 
eral description of a large class of social systems object of the modeling approach 
and outlines the modeling strategy which generates the mathematical structure 
we look for. Subsection 3.2 defines the representation of the system by link¬ 
ing at each functional subsystem a probability distribution over the microscopic 
state. In subsection 3.3 the aforementioned structure is derived. Finally, some 
reasonings on the modeling of interactions are discussed. 

3.1 An heuristic description of social systems 

Let us now define precisely the class of systems which is object of the modeling 
approach. In detail, we consider a network of a number m of interacting nodes, 
where in each node individual entities are grouped into n functional subsystems. 

According to m, the hallmarks of the modeling approach can be summa¬ 
rized as follows: 

• The entities which comprises the system are referred to as active particles 
and are assumed to be distributed in a network of nodes labeled by the 
subscript i = 1,... ,m. Therefore, the role of space variable is confined 


to the network, namely a continuous variation in space is replaced by 
migration dynamics across the nodes. 

• Active particles are partitioned into functional subsystems, where they 
collectively develop a common strategy. These are labeled by the subscript 
j = 1,... ,n. Therefore the total number of subsystems is N = mn. 

• The strategy is heterogeneously distributed among active particles and 
corresponds to an individual state, defined activity. Hence the state of 
each functional subsystem is defined by a probability distribution over the 
said activity variable. 

• Active particles interact within the same functional subsystem as well as 
with particles of other subsystems, and are acted upon the external ac¬ 
tions. These interactions are generally nonlinearly additive and are mod¬ 
eled as stochastic games, meaning that the outcome of a single interaction 
event can be known only in probability. 

• The evolution of the probability distribution is obtained by a balance of 
particles within elementary volumes of the space of microscopic states, 
the inflow and outflow of particles being related to the aforementioned 
interactions. 

These hallmarks are quite general and need to be more precisely defined de¬ 
pending on the system under consideration. In particular, the following aspects 
have to be accounted for: 

i) Definition of the functional subsystems and modelling of their dynamics across 
nodes; 

ii) Modelling of the transitions across functional subsystems; 

iii) Presence of leaders, which should be grouped in specific functional subsys¬ 
tems; 

iv) Characterization of the so-called external actions and of the dynamics in¬ 
duced by them; 

Remark 3.1 A vast literature has been recently developed on the theory of net¬ 
works and applications, for example [II3 [22f ■ 

This present paper refers only tangentially to such a theory, while some theo¬ 
retical tools of evolutive game theory are used to model interactions also across 
the network. 

3.2 Functional snbsystems and their stochastic represen¬ 
tation 

The overall state of the internal system is delivered by a probability distribution 
over the activity variable in each node and for each functional subsystem, here¬ 
inafter denoted by the acronym ij-FS, while the generic one is simply denoted 
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by FS. In detail, the probability distribution reads 

: [0,T] X Dii —>• K_|_, (1) 

such that moments can provide a description of quantities of the system at the 
macroscopic scale. In detail the zeroth order moment gives the number density: 

riij [fij] (t) = E°. [fij] (t) = [ fij {t, u) du, (2) 

Jdu 

which gives the number of active particles in each ij-FS. Higher order moments, 
corresponding to additional macroscopic variables, can be computed as follows; 

lEy(i)[/zj] = —T^-rTTY / uPfij{t,u)du, p = l,2,..., (3) 

where the computing of macroscopic quantities is meaningful only if the number 
of active particles in the FS is a strictly positive defined quantity. 

The activity is represented by scalar variable whose domain of definition 
Du is typically [—1,1] or even the whole real axis for probability distributions 
which decay to zero at infinity. Positive values indicate an expression coherent 
with the activity variable, while negative values expressions in opposition. The 
interval [—1,1] is used whenever a maximum admissible value of the activity 
variable can be identified so that u is normalized to such value. 

More in general the activity variable can be a vector, but we remain here 
within the relatively simpler case of a scalar variable, which makes light the 
formalism of the mathematical approach. 

The actions exerted by the external system are supposed to act at the same 
scale of the system and over the same domain of the activity variable by agent 
systems. These actions are assumed to be known functions of time and of the 
activity variable defined by the following structure: 

(Pijit,u) = Sij{t)ipij{u) : [0,r] X Du ^ K+, (4) 

where models the intensity of the action in the time interval [0, T], while 'ipij 
models the intensity of the action over the activity variable. 

3.3 Mathematical structures and models of interactions 

The mathematical theory of active particles m leads to the derivation of math¬ 
ematical structures, which describe the evolution in time of the probability dis¬ 
tribution over the activity variable for each functional subsystems. The result, 
which consists in a system of n x m differential-integral equations, is obtained 
by a balance of particles within the elementary volume [u, m -I- du] of the space 
of the microscopic states. This approach is documented in the already cited 
survey m as well as in various papers concerning specific applications such as 
the immune competition |341 166] or learning dynamics [56] . Therefore, only a 
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concise description of the technical aspects of the approach will be given here 
referring to the existing literature for further details. 

The dynamics of interactions can be obtained by assigning to each particles 
of the internal ij-FS the following distinguished roles: Test particle of the ith 
node and jth FS with activity u, at time t; Candidate particle of the ith node 
and jth FS with activity it*, at time t, which can gain the test state u as a 
consequence of the interactions; Field particle with activity it*, at time t, which 
triggers the interactions of the candidate particles. Hereinafter, the term (ij)- 
particle is used to denote a particle in the ijth FS. It is assumed that the overall 
state of the system is delivered by the probability distribution linked to the test 
particle. 

Moreover, the following interactions occurring at different scales are consid¬ 
ered: 

• Microscopic scale interactions: Individual based interactions between 
particles belonging to the same FS or to different FSs. 

• Microscopic-macroscopic scale interactions: Interactions between 
particles and FSs viewed as a whole being represented by their mean value. 

The mathematical framework describing different kind of interactions will 
be specified, for each type of interactions, by means of two terms: 

• The interaction rate which describes the rate of interactions involving 
particles; 

• The transition probability density which describes the probability density 
that a candidate particle falls into the state of the test particle after the 
interaction. 

Details on the modeling of these interaction terms will be given in the next 
subsection. Supposing now that these terms can be defined by appropriate 
models, the following balance of number of particles can be used to derive the 
mathematical structure: 

Dynamics of the number of active particles 

= Net flux due to ms-interactions between active particles 
+ ms-interactions between active particles and external actions 
-h mMs-interactions between particles and functional subsystems, 

where the net flux denotes the difference between inlet and outlet fluxes, while 
“ms” stands for microscopic scale and “mMs” stands for microscopic-macroscopic 
scale. 

Detailed calculations are not reported here as the interested reader can find 
them in some recent papers such as [40l [66] and [3l] . The result of the said 
calculations is the following structure: 

dtfij{t,u) = Aij[i]{t,u) + B^j[i,ip\{t,u) + Cij[i]{t,u), (5) 
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where 


m n p p 
„ I- 1- 7- . J J D„ xD„ 


p^h—1 q,k—l ' 

X fpq{t,U^) fhk{t,u*)du^du* 


m n p 

i ’ 7 *'*/[f](u, u*) fhk{t, U*) du* 

h=l fc=l 


( 6 ) 


Bij[i]{t,u) = XX// ^i.pq[^,'^]{u^,V*)Bpq{ij)[^,if\{u^, ^ u\u^,v*) 

p=l q=l d JDuXDu 
X fpq {t, U*) (Ppq{t, U*) dUt. du* 

- fijit^u) f lJ.tj[{,ip]iu,v*)ipij{t,u*)du*, (7) 

d Du 

m n p 

X X / ^P9 [fKw*)C“(7j)[f](M, M|u*,E7ifc[f])/p9(t,U*) 

[f] 

m n 

/.(t,.)XX 4'=[f](M)E^fe. (8) 

/l —1 

Moreover, the interaction rates, which model the frequency of interactions for 
each type of interaction have been denoted by 77 ^^, fipq and . The transition 
probability densities, which denote the probability density that a candidate pq- 
particle shifts, after an interaction with a field hfc-particle, to the state of the test 
ij-particle have been denoted by Apq, Bpq, and . The output of interactions 
has been put into round brackets. Square brackets have been used to denote 
dependence on the distribution functions (nonlinear interactions). 

Remark 3.2 It has been assumed that each node includes all FSs. This as¬ 
sumption is formal if the initial status of some FSs is equal to zero and remains 
like that for t > 0. On the other hand, such status can take values for t > 0 due 
to migration phenomena. 

Remark 3.3 This structure includes some specific systems. As an example, 
the case of a closed system can be considered by setting ip = 0. The structure 
simplifies as follows: 

dtfij(t,u) = +Cij[i]{t,u), (9) 

On the other hand, a possible alternative to the structure can be obtained by 
modeling an external action Kij(t,u) acting over each functional subsystem, 
which yields: 

dtfij{t,u) + du[Kij{t,u)fij{t,u)\ = + Cij[f]{t,u). (10) 


C^j[i]{t,u) = 

X 
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Or hy supposing that the external actions acts over the system by the mean value 
of ipij so that the term Bij is modified aecordingly. Further generalizations can 
be achieved by taking into account also macroscale actions by each node as a 
whole and by the whole network. 

Remark 3.4 The output of interactions between two interacting entities can be 
conditioned by their states at the microscopic scale and by the probability distri¬ 
bution over these states. As previously mentioned, interactions are linear if the 
output is eonditioned only by the microscopic states, while are nonlinear when 
the output is conditioned by the probability distributions or by their moments. 
The various sources of nonlinearity featured by the eguation above are critically 
analyzed in J33\/ . 

Remark 3.5 The term stochastic games is used since the state of the interact¬ 
ing entities, as well as the output of their interactions, are known in probability. 
Moreover, the term evolutionary is used to account for the dynamics in time of 
the rules that drive interactions. Since these games are casted into a differential 
system, we will talk about stochastic evolutionary differential games. 

Remark 3.6 Several applications, as examples suggest to use discrete 

variables at the microscopic scale. In fact, in some specific cases the state of the 
active particles is more precisely identified by means of ranges of values rather 
than by a continuous variable. Technically, this means that the distribution 
functions, for each FS, have to be regarded as discrete distribution functions, 
while the interaction terms map discrete variables rather than continuous ones. 
The mathematical structures are readily obtained by specializing the integrals 
over the microscopic states as sums over discrete states. We skip over the tech¬ 
nical, it might be tedious, derivation of a general expression of a mathematical 
framework for systems with discrete states, as it is more practical focusing such 
a problem on specific models rather than on a general case. This topic will be 
treated, at a practical level, later in Subsection 5.2 when the general approach 
will be focused on a case study. 

3.4 Some reasonings towards modeling interactions 

Specific models of social dynamics can be obtained by inserting models suitable 
to depict interactions at the microscopic scale into the structure ([5]). Various 
recent papers have contributed to this topic. For instance [32l |36l |69j have 
shown how interactions can be modeled by games where the output of the 
interactions is conditioned not only by the state of the interacting entities, 
but also by the probability distribution over such states. Moreover, models 
of multiscale micro-macro interactions have been proposed in [100] concerning 
the modeling of migration phenomena, in [101] on opinion formation in small 
networks, and in [32] focused on criminality dynamics. These multiscale features 
deserve attention in several fields of social dynamics. 

The aforementioned literature suggests to critically revised the whole sub¬ 
ject. Such a revisiting is needed to derive tools that are valid for a large variety 
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of social systems, as well as to deal with the actions externally exerted on the 
system, which is a topic only tangentially treated in the literature although very 
important in the modeling of social systems. 

Bearing all above in mind, let us consider, separately, the modeling of the 
interaction rates and of the transition probability density. 

3.4.1 Interaction rate 

The modeling of interactions requires the definition of different concepts of dis¬ 
tance between interacting entities [32) . In details: 

• Microstate-microstate distance refers to the states at the microscopic scale 
between candidate and field particles. This distance can be |u* — u*\. 

• Microstate-macrostate distance refers to the interaction of a candidate 
particle with state u* with the group of particles belonging to a FS with 
mean value of the activity E. This distance is given by |u* — E|. 

• Affinity distance refers to the interaction between active particles char¬ 
acterized by different distribution functions. This distance is introduced 
according to the general idea that two systems with close distributions are 
affine and is given by \\fpq — fhk\\, where || ■ || is a suitable norm to be 
chosen depending on the physics of the system under consideration. 

• Hierarchy distance whenever a hierarchy exists among groups of interest 
to distinguish the more and less relevant ones. 

The overall distance, let us call it social metrics, is a weighted sum, which 
replaces the geometrical metrics, of all distances. In general, the interaction 
decays with the distance, where heuristic assumptions lead to a decay described 
by exponential terms or rational fractions. 

3.4.2 Transition probability density 

The dynamics of interactions can be modeled by theoretical tools from evo¬ 
lutionary and behavioral game theory (57) [HHl HD 1112111191112011121) . which 
provides features to be introduced into the general mathematical structure in 
order to obtain specific models. Additional tools are given by learning the¬ 
ory and evolutionary games m, within the broader context of statistical 
dynamics and probability theory. For instance, the following types of games are 
often used: 

• Competition (dissent): The interacting particle with higher status in¬ 
creases its status by taking advantage of the other with lower status. 
Therefore the competition is advantageous for only one of the two players 
involved in the game. 


14 












• Cooperation (consensus): The interacting particle shows a trend to 
share their microscopic states by decreasing of the difference between the 
interacting particles’ states, due to a sort of attraction effect P1ITM|. 

• Learning: One of the two particles modifies, independently of the other, 
its microscopic state, while the other reduces the distance by a learning 
process. 

• Hiding/chasing: One of the two particles attempts to increase the dis¬ 
tance from the state of the other one [hiding)^ which conversely tries to 
reduce it (chasing) [5^ . 

Remark 3.7 All aforesaid types of games can occur simultaneously in a general 
eontext of heterogeneous particles. In some cases, see the oecurrenee of 
one of them is ruled by a threshold on the distance between the states of the 
interacting particles. Such a threshold is, in the simplest case, a constant value. 
However, it is worth mentioning that a recent paper \69f has shown that the 
threshold can depend on the state of the system as a whole and that can have 
an important influence on the overall dynamics. Thus dynamics refers to the 
interaction of more than one type of dynamics as treated in m from the view 
point of social economy, where political choices are related to economical growth. 
Further applications, in addition to the already cited paper are proposed 
in iSblj and 

Let us now consider the above remark referring more in detail [69) . where a 
parameter has been introduced to account for the “selfishness” of the society. 
Such parameter is related to the evolution equation describing the dynamics of 
the aforementioned threshold. A detailed study of the model shows that high 
values of the selfishness in the wealth distribution ends up with an overall re¬ 
duction of the total wealth, which is shared only by small groups, while the 
great part of the population is confined in poverty. The mathematical structure 
has been modified in [69j by linking the equation for the dynamics of the dis¬ 
tribution function to an equation for the threshold triggering interactions. The 
formal structure in the general case is as follows: 


fij (t, u) — [f, /r] (t, n) -f Cij [f, /x] (t, u ), 




( 11 ) 


where p. denotes the set of all thresholds /x^ . A deep insight into the role of 
selfishness and the complex problem of its “measurement” is proposed in |124j . 
The effect of network topology on wealth distributions is studied in m- 

Based on simulations. Paper [53] predicts that in an excessively selfish soci¬ 
ety wealthy classes initially enrich themselves but eventually the whole society 
moves toward poverty since the middle classes, which generally are the produc¬ 
ers of positive growth of the market, are strongly weakened. 
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Remark 3.8 The interaction domain Du has been supposed, for simplicity, 
constant. However, it can he assumed to be depending on the distribution func¬ 
tion. A possible example is that of i38\/ based on the conjecture that sensitivity is 
related to a critical size rather than to the whole perception domain. This concept 
deserves attention to be related to specific models. This remark accounts for \38\/ . 
where it is suggested that interactions occur in a domain a sensitivity domain 
^ C Du, corresponding to a critical number Uc of field particles. Integration of 
fi over the activity variable in a domain fl = [u —Sm[f',nc],u-\-SM[f',nc]], which 
can be called the topological domain of interaction, can lead to compute Sm and 
sm, where Sm,SM > 0. However, the solution is unique only in some special 
cases. For instance, when u is a scalar defined over the whole real axis, and the 
sensitivity is symmetric with respect to u. On the other hand, if u is defined in 
a bounded domain or the sensitivity is not symmetric, additional assumptions 
are required. 


4 A Critical Overview of the Existing Literature 

The mathematical structure proposed in the preceding section has been al¬ 
ready used, simplified to special cases, by some models known in the literature. 
Therefore, it is worth developing an overview of these particular models to un¬ 
derstand which type of innovation can be achieved by a more general approach. 
This analysis can contribute to show how the various approaches are related to 
the complexity features of social systems. These two topics are presented in the 
next two subsections as an answer to the third key question: 

How far the existing literature has exploited a “mathematical 

structure ”? 

Before presenting a survey of the existing literature, let us mention that the 
search of a structure, which is here viewed as a necessary step of the model¬ 
ing approach, can have a broader meaning in various fields of life sciences as 
observed by Gromov [S3, who stresses that a deep analysis of mathematical 
structures might offer an overview much reacher than that foreseen by the ini¬ 
tial application. Referring to the specific structure under consideration, suitable 
developments such as the introduction of proliferative and destructive interac¬ 
tions have generated models in biology |40j . while the introduction of a space 
structure have lead to models of vehicular traffic and crowds [SllSilSS]- 

4.1 An overview of models 

The application of kinetic type equations with interactions modeled by theo¬ 
retical tools of game theory was introduced in m and subsequently paper [ll] 
generalized such a structure to model large social systems with discrete states 
at the microscopic scale. The authors used games with linear binary interac¬ 
tions and derived a system of ordinary differential equations suitable to describe 
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the time dynamics of a discrete probability distribution corresponding to the 
discrete micro-states. Actually, the authors used a simplified version of a frame¬ 
work proposed in m to model complex multicellular systems. The application 
developed in [33] refers to the dynamics of wealth distribution and is based 
on consensus/dissent games triggered by a constant threshold of the distance 
between the states of the interaction pairs. 

The same approach was subsequently extended with minor technical mod¬ 
ifications to various types of applications by the authors and their coworkers. 
Without claim of completeness, examples include opinion formation [45] . tax¬ 
ation dynamics with some social implications |3S], value estimates [5S]. Vari¬ 
ous innovative developments followed these pioneer papers. The main sequen¬ 
tial steps of the aforementioned developments are summarized in the following 
looking ahead to new research approaches, which can be achieved by use of the 
structure of Eq. 

Opinion formation is an important ground of application of mathematics to 
the modeling of social sciences, hence several valuable papers have been devoted 
to this topic, for example fni 1301 mHIM IMl [HT] have developed interesting 
models based on kinetic theory methods. These methods have also been applied 
to modeling financial markets Esioa. 

A pioneer approach to systems theory: Various authors have proposed 
different approaches to a systems theory of social dynamics, for instance [UliOj. 
However, the first contribution along the conceptual line of this present paper 
has been given in |5], where the concept of scaling and functional subsystems 
have been proposed referring to a variety of social systems as well as the obser¬ 
vation that modeling the interactions of different types of dynamics is needed. 
This is, indeed, a topic of great interest in socio-political sciences [47] . 

On the introduction of nonlinear microscopic interactions: Nonlin¬ 
earity for interactions at the micro-scale have been introduced in [36] and [69] . 
Both papers analyze the interplay of different dynamics, namely wealth dis¬ 
tribution versus support or opposition to Governments, or selfishness in wealth 
dynamics as a source of overall wealth dissipation in nations. The role of nonlin¬ 
earity of interactions was clearly put in evidence in [33] . An interesting reference 
for social aspects of selfishness and wealth distribution is given by |116j . based 
on previous studies in the field [85] 110311104] . 

Dynamics over small networks: The modeling of this type of dynamics 
bv kinetic theorv methods was introduced in [100] and m referred to specific 
applications, such as migration phenomena. These two papers have introduced 
a new type of nonlinearity under the assumption that individuals are sensitive to 
the mean value of each functional subsystem, which depends on the distribution 
function. Indeed, this is the two scales micro-macro interaction presented in our 
paper. An alternative approach to modeling migration dynamics are presented 
in m and [13]. The concept of social distance, introduced in Subsection 4.3.1, 
is one of the main tools toward the study of networks. Such a distance can 
be viewed as a stochastic quantity as it depends on a probability distribution. 
Therefore the structure of a social network evolves in time and is not rapidly 
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related to the physical localization. 

Mutations and selection: A post Darwinian dynamics consisting in mu¬ 
tations followed by selection plays an important role in biology, specifically 
in the immune competition |341 1401 I66j . An analogous dynamics appears in 
socio-economical systems, where new groups of interest can be generated in the 
dynamics, for instance by aggregation of different groups, and subsequently they 
might either expand or disappear in a competition somewhat mediated by the 
environment where these groups act. These concepts have been introduced in 
some recent papers, for instance |64j in a behavioral theory of urbanism. 

Control problems: An important field of investigation is the development 
of a control theory of social systems. This topic has been developed in the 
study of swarm dynamics |551159j and recently focused more precisely on the 
dynamics of social systems |58| . The mathematical structures proposed in this 
paper include the presence of external actions to be properly studied to drive 
the dynamics toward optimal behaviors. 

4.2 Descriptive ability of the mathematical structures 

A critical analysis of the capability of the mathematical structure ([5]) to capture 
the complexity features presented in Section 2 is a necessary step toward the 
systems approach. In detail: 

• Ability to express a strategy: The microscopic state includes the 
activity variable, which is deemed to model the strategy expressed by each 
functional subsystem. 

• Heterogeneity: The heterogeneous behavior is accounted for by the de¬ 
scription of the system by means of a probability distribution linked to each 
functional subsystem. 

• Nonlinear interactions: Nonlinearity of interactions has been already 
discussed in the preceding section. As a result linear additivity cannot 
be applied. Moreover, we observe that an additional source of nonlinear¬ 
ity is related to non-locality of interactions between functional subsystems 
localized in different nodes of the network. 

• Learning and adaptation: The ability of living systems of learning from 
past experience can be mimic by modeling interactions based on rules that 
evolve in time and include a continuous adaptation to the changing-in-time 
environmental conditions. 

• Selection and evolution: The onset of functional subsystems is mod¬ 
eled by the transition probability density, where interactions include the 
formation of groups of interest, which in turn can generate new groups 
more suited to an evolving social and economical environment. Selection 
is modeled by interactions with the outer environment. 
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It is worth noticing that the descriptive ability of the complexity features 
is only potential and it should be practically implemented in the derivation of 
models when theoretical tools of game theory are used to describe interactions. 
Although models rapidly reviewed in Subsection 4.1, have been successfully 
depicted some collective behaviors observed in reality, it can be remarked that 
one of the limits of the existing literature is that it deals with closed systems, 
while interactions with the outer environment has not been taken into account. 
Moreover, the role of networks in determining the overall dynamics has been 
treated only tangentially. Therefore, a great deal of work is still waiting for being 
developed. Hopefully, the next sections will contribute to link the conceptual 
ideas proposed until now to real applications. 


5 Toward a Systems Approach to Behavioral So¬ 
cial Dynamics 

This section shows how the mathematical tools presented in Section 3 can be 
applied to modeling and simulation of social systems. The aim consists in 
providing an answer to the fourth key question: 

Can a general systems theory be developed and how can it 
contribute to the interpretation of the so-called “big-data”? 

With this aim in mind the contents of this section can be presented. Sub¬ 
section 5.1 provides the general hallmarks of the theoretical approach. Then, 
Subsection 5.2 provides a critical overview on Monte Carlo computational meth¬ 
ods, which appear to be the most efficient method for the class of equations 
proposed in this paper. Finally, Subsection 5.3 shows how the system approach 
can contribute to the use and interpretation of repositories of the so-called big 
data toward the interpretation of the dynamics of social systems. This section 
will remain at a theoretical level, while a specific case study is presented in 
Section 7, where various concepts are referred to a well defined model. 

5.1 Hallmarks of a systems theory 

Some guidelines towards a systems approach have already been anticipated in 
Subsection 3.1, however not yet referred to any mathematical formalization. 
The mathematical tools derived in Section 4 allow us to proceed in this crucial 
step. Therefore, a revisiting of the previous concepts is proposed as follows: 

1. The domain of the space variable where the system is located is subdivided 
into interconnected sub-domains called nodes, while the set of all nodes 
and their connection is called network. 

2. The individual entities, called aetive particles, are aggregated into different 
groups of interest called funetional subsystems. 
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3. The active particles aggregated in a functional subsystem express a com¬ 
mon strategy (or expression of interest) called activity. 

4. The heterogeneity of the particles’ activity within each functional sub¬ 
system is expressed through a probability distribution. The overall set of 
probability distributions represents the dependent variable of the dynam¬ 
ics. 

5. Modeling is required for interactions between active particles within the 
same functional subsystem and with particles of other subsystems, as well 
as for external actions. All these interactions can induce transitions across 
functional subsystems. 

6. The equations modeling the dynamics of the probability distributions over 
the micro-states are obtained by a balance of particles within elementary 
volumes of the space of microscopic states, the inflow and outflow of par¬ 
ticles being related to the aforementioned interactions. 

7. These equations can be solved by suitable computational methods to ob¬ 
tain both the probability distributions and the moments which provide an 
information of interest on macroscopic states. 

Remark 5.1 The following speeific milestones of the approach: Subdivision, 
activity variables, external actions, network, etc. are not constant features of 
the overall system. In fact, these variables depend on the specific social dy¬ 
namics object of study. Therefore, different investigations would lead different 
characterization of the overall system and hence of mathematieal models. 

5.2 Computational methods 

The development of the numerical methods for solving the class of equations 
proposed in this paper should take into account that the dependent variables 
are probability distributions over the activity variable depending on time. De¬ 
terministic solutions can be achieved by approximating the said distribution at 
a number of collocation points over the activity variable. This transforms the 
original integro-differential system into a coupled system of ordinary differential 
equations corresponding to the discrete values of the probability distributions in 
each point of the collocation. Classical numerical techniques for solving ordinary 
differential equations can then be employed. Examples of this method of solu¬ 
tion applied to equations whose mathematical structure closely resembles the 
one presented in the present paper can be found in [711H1IH3]- Notice should 
be made that the resulting system of ordinary differential equation can be usu¬ 
ally very large for a network where in each node various functional subsystems 
interact. Indeed this is the case of the developments on the modeling of the 
criminality dynamics reported in the following section. Although parallel com¬ 
puting can be used to alleviate the computational burden [zaize], alternative 
numerical approaches are certainly of interest. 
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According to the authors’ bias, Monte Carlo simulation [5711115) methods are 
ideally suited for solving the class of equations under consideration. The basic 
idea consists in representing the distribution function by a number of compu¬ 
tational particles. These particles move through the network, unless migration 
dynamics is disregarded, and interact according to stochastic rules derived from 
the kinetic equation, Eqs. ([S]). The macroscopic fields are obtained through 
weighted averages of the particle properties. Compared to deterministic meth¬ 
ods of solution, the application of stochastic methods provides some important 
advantages such as the computational efficiency and the possibilities to easily 
account for sophisticated individual decision processes. 

It is worth noticing that a distinction is sometimes made between simulation 
and Monte Carlo, being the former the direct coding of a natural stochastic pro¬ 
cess while the latter the solution by probabilistic methods of non probabilistic 
problems. However, often such a distinction cannot be maintained. A relevant 
example is provided by the Direct Simulation Monte Carlo (DMSC), which is 
by far the most popular simulation method for classical kinetic theory applica¬ 
tions. Although it has been introduced based on physical reasoning [35], DSMC 
has been later proved to converge, in a suitable limit, to the solution of the 
Boltzmann equation [133] . Therefore the same computer code can be regarded 
simultaneously as a direct physical transcription of the dynamics of a rarefied 
gas and as a stochastic solution of the Boltzmann equation. 

5.3 On the use and abuse of the so-called big data 

This subsection aims at showing how the systems approach can contribute to 
the use of the so-called big data toward predictive purposes on an ongoing 
observed dynamics. As it is known, thanks to the advent of information and 
communication technology, as well as to the increasing amount of (big) data 
is collected everyday. These data contain information about the behavior of 
individuals (such as their movements, their networks of friends, the sharing 
of opinions via social networks, etc) and are increasingly becoming a valuable 
source of information to validate models around the behaviors of individuals. In 
general, the term big data is used to define data sets so large or complex that 
traditional data processing applications are inadequate. An important reference 
is given by the report m by paper [68] . while additional information can be 
obtained by the reports posted in the WEB sites [55] . 

The problem consists in understanding how these data can be used for pre¬ 
dictive purposes or for model validation. Currently, the various methods to 
treat these large amounts of data are being developed define an emerging - so 
called - data science which aims at improving the decision process toward cost 
reductions and reduced risk. Rather than treating the technical problem of data 
compression, we suggest, as a possible perspective, the use of big data within the 
systems approach proposed in this paper focused on the analysis of an observed 
dynamics where the main problem consists in selecting specific actions to drive 
the dynamics toward an optimal trend. This target refers to a large number of 
data collected in a repository database. The approach can be developed along 
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the following sequence of actions: 

1. Classification of the data according to the specific features of the systems 
approach, namely by functional subsystems and activity variables. These 
features are used to compare different dynamics independently on external 
actions. 

2. Define an appropriate metrics suitable to compare the specific features 
of the system under consideration with similar dynamics stored in the 
aforementioned database. Such a metrics defines the distance between 
the selected dynamics Sd and the stored one Ss- Us'd — S'r||, where, after 
appropriate normalization, ||Sd —Sr.|| =0 means maximal similarity, while 
||Ss — Sr-ll = 1 means maximal distance. 

3. Select in the big data repository the closest Ss to the system under con¬ 
sideration and provide for each of them a score V € [0,1], where D = 1 is 
the best score, while V = 0 means the worst one. 

4. Select in the database the optimal dynamics, namely the optimal external 
actions, by using a convex combination of the distance and the score. For 
instance: 

J = a(l- ||Sd-S,||) + (l-a)D, 
where the optimal selection corresponds to the highest value of JV. 

5. Apply to the real systems the external actions corresponding to the dy¬ 
namics in the database corresponding to the system selected according to 
item (4). 

Remark 5.2 The milestones (l)-(5) allow to apply a decision process related to 
the selection of actions appropriate to control the dynamics of systems by selec¬ 
tion actions appropriate to control an ongoing dynamics. The actual expression 
of of the objective functional can be further revised by enriching it of additional 
variables according to the specific system which is treated. The method is general 
and can go beyond applications in social dynamics. 

6 Reasonings on a Case Study 

Let us now show how the general systems approach proposed in the preceding 
section can be applied in a specific case study. In detail, we consider the onset 
and development of criminality in a society, where this dynamics is contrasted 
by security forces, such as intelligence and police. 

The choice of the case study refers to a mathematical literature, which is 
constantly growing due also to the impact that this topic has on the wellbeing 
and security of citizens. Useful references are given by the following essays [4^ 
m [71 [HHi [m, which are brought to the attention of the reader as examples 
without aim of completeness. 
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Different mathematical approaches have been used corresponding to different 
representation scale and tools to model interactions. Nonlinear dynamical sys¬ 
tems have been develooed in |113| . while the use of comnartmental models m 
allows modeling transitions across what we have called functional subsystems, 
and predicting, at least in some specific cases, tipping points. The approach 
by partial differential equations, which include not only number of individuals, 
but also their localization has been proposed in various papers, for example [43] . 
where accounting for the localization includes the modeling of pattern formation 
of criminality and its propagation in the territory. The dynamics of shocks and 
acceleration waves are studied in [128] . 

This section illustrates how the hallmarks of the system theory previously 
discussed can be applied to the model presented in [32], which shows how the 
contrast actions are related on one side to the social state of citizens, such as 
wealth and culture, and, on the other hand, to the ability of security forces 
based also on their training. More in general, the interplay between wealth 
distribution and unethical behaviors is object of growing interest of researchers, 
who operate in the field of social sciences [104] , while it is becoming a new field 
of investigation of applied mathematicians [36l |69] . 

The contents are presented in three subsections. Subsection 6.1 briefly sum¬ 
marizes the model proposed in [32] . Subsection 6.2 shows how a model with 
discrete states can be designed in agreement with Remark 3.6. Finally Subsec¬ 
tion 6.3 proposes some possible developments. 

6.1 On a model of social criminality 

The mathematical model studied in [32] accounts for the following interaction 
dynamics: Susceptibility of citizens to become criminals, susceptibility of crimi¬ 
nals to reach back the state of normal citizens, learning dynamics among crimi¬ 
nals, motivation/efficacy of security forces to catch criminals, learning dynamics 
among security and interplay with wealth dynamics. 

In more details, the following features of the model proposed in [32], are 
briefly analyzed: 

Functional subsystems and activity: The following subsystems are con¬ 
sidered: Citizens, criminals and security forces, while the activity variables ex¬ 
pressed by them are, respectively, the following: Wealth, criminal ability and 
detective ability. 

Mathematical structure: The mathematical structure used in [32] is the 
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following: 


dtfi{t,u) = Ji[{]{t,u) = 


= ^ / / r]hk{u*,u*)Bl,,{u^. ^ u\u^,,u*)fh{t,u^)fk{t,u*)du^.du* 

h,k=i 

3 

/ r]ik{u,u*) fk{t,u*)du* 

k=l 

J Di 


-^i(u,Ei)/i(f,u), 


where : [0,T] x Di ^ K.+ , i = 1,2,3, while r]hkiu*,u*) and /ii(tt*,Ei) are 
the encounter rates between a candidate h-particle and a field /c-particle and 
between a candidate z-particle and the mean activity, respectively. Moreover, 
—>■ and —)• M|u*,E;i) are the probability density for the 

state transition of individual due to interactions between particles and between 
particle and the mean activity, respectively. Such a structure is a particular case 
of the general one reported in Eq.(3.5)-(3.9). It is worth noticing that external 
actions are not accounted for and space dynamics is disregarded. Therefore this 
model refers to the simple case of a closed system which comprises only one 
node. 


Modeling encounter rates: Citizens with closer social states interact more 
frequently; unexperienced lawbreakers tend to expose themselves more frequently; 
experienced detectives are more likely to chase less experienced criminals. 

Modeling individual based interactions: Citizens are susceptible to be¬ 
come criminals, motivated by their wealth state, thus they can mutate into a 
new functional subsystem. Criminals interact among themselves resulting in 
a dynamics by which less experienced criminals mimic the more experienced 
ones. Criminals chased by detectives are constrained to step back decreasing 
their activity value as the price to be paid for being caught. At the same time, 
detectives gain experience from a well-done job increasing their activity. Due to 
this action, criminals are induced to return to the state of normal citizens with 
probability which increases with decreasing values of their level of criminality 
and increasing values of skill of detectives. 

Interactions between individuals and the whole functional subsystem: 

Both criminals and detectives interact with the mean value through the mean- 
micro state distance within their own functional subsystem, in both cases the 
interaction rate increases with the distance between the individual state and 
the mean value. The dynamics is such that only those who are less experienced 
than the mean tend to learn and move towards it. Likewise, detectives show a 
trend toward the mean value. 

For additional technical details on the terms modeling interactions, namely 
interaction rates and transition probability densities, we refer the reader to |32j . 
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where a number of simulations are carried out which indicate that an equal dis¬ 
tribution of wealth leads to a slow growth, or even a decrease of criminals, while 
the opposite trend is observed in the case of unequal distributions. Concerning 
the role of detectives the model shows that the expertise which they acquire due 
to training more than their number has a positive effect on the control of the 
number of criminals and their aggressive behaviors. 


6.2 On a model with discrete states 

As pointed out in the Remark 3.6, the use of models with discrete states is moti¬ 
vated by the difficulty of dealing with a continuous activity variable. Moreover, 
in several cases, the discrete approach can simplify the modeling of interactions 
at the microscopic scale and even enriches the predictive ability of the model. 
This subsection aims at showing how a model with discrete states can be de¬ 
rived under the same phenomenological assumptions proposed in the preceding 
subsection. 

As in Section 6.1, the functional subsystems and related activity variables are 
the following: i = 1 : citizens and wealth; i = 2 : criminals and criminal ability; 
i = S : detectives and their skill. However, we here introduce a discrete activity 
variable, labeled by the subscript cr = 1,..., m with ui = 0 and Um = 1- The 
abbreviation i; cr-particle is used for a particle belonging to the f-th functional 
subsystem with activity u^. The representation of the overall system is delivered 
by the discrete probability /f = /f(t) : [0,T] ^ ]R_|_, while weighted sums 
replace the weighted integrals to obtain variables at the macroscopic scale. 

The actual derivation needs defining the interaction terms in the new context 
as well as a new mathematical structure. In detail: 


• is the encounter rate between a candidate h; p-particle and a field 
k; g-particle. 

• is the encounter rate between a candidate h;p-particle and the mean 
activity within its functional subsystem. 

• cr) is the discrete transition probability that a candidate h; p-particle 
ends up into the state of the test i; cr-particle due the interaction with a 
field k; g-particle, such that 

• <j){up, E^) is the transition probability that a candidate /i;p-particle 
ends up into the state i; cr-particle due the interaction with the mean 
activity value Eh- such that that 


The transition probabilities satisfy, for all type of inputs, the following nor¬ 
malization 


3 m 


m 


E E ^) = E = 1. (12) 
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Derivation of a mathematical structure The balance of particles in the 
elementary volume of the space of micro-states leads to the following structure: 


dtf[it) = jnm = 


3 m 3 m 



h,k—l p,q—l 


k—1 q—1 


3 m m 


m 




(13) 


h=lq=l 


The derivation of the mathematical model is obtained within the general 
structure given by Eo. lfT^ by particularizing the terms 77, p,, B and At. 

Encounter rate Let us consider the encounter rates and The modeling 
approach is based on heuristic assumptions that let us quantify the frequency 
of interactions, depending on the micro-states and distribution functions of the 
interacting particles. The result, where candidate particles are represented by 
a square, while field particles are represented by a circle, is the following: 

• [^ O @: Closer social states interact more frequently r]ii = {1 — \up — Uq\) 

• O @: Experienced lawbreakers are more expected to expose them¬ 
selves 7722 = ff(Up + Uq)\ 

• @: Experienced detectives are more likely to less experienced crim¬ 
inals hunt 772! = ((1 “ ^p) + '“?) i 

• O E2: Criminals interact with the mean value through the mean-micro 
state distance 772(17*,E2) = 77°|u* — E2I; 

• [3] O E3: Detectives interact with the mean value through the mean-micro 
state distance 773(17*,E3) = 77°|i7* — E3I. 

Transition probability: The social structure of the population i = 1 is as¬ 
sumed to be fixed, namely, the time interval is sufficiently short, that the wealth 
distribution is constant in time. 

• The probability that citizens become a criminals, increases with their de¬ 
creasing wealth state referred to that of other citizens: 


= OA (1 - Up)Uq, Bll{\\) = (1 - OA (1 - Up)Uq). 


• Less experienced criminals mimic the more experienced ones, moreover 
also interaction with less experienced lawbreakers increases the level of 
criminality: 


^22(2(P + 1)) = /3(1 - Up)Uq, Bll{2p) = (1 - ^(1 - Up)Uq). 
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• Criminals are constrained by detectives to step back decreasing their ac¬ 
tivity value as the price to be paid for being caught: 

p = 1 : 1) = asim- l)up, Sf|(2,1) = (1 - asim - 1)^^); 

P > 1 : ^23(2,^- 1) = aBirn-p)up, i 3 f|( 2 ,p) = (1 - aBim - l)up). 

• Detectives gain experience from their activity and increase their skill: 

^32(3,P + 1 ) = 7^5(1 - Up), (1 - 7 Mq(l - Up)). 

• The level criminality and the skill of detectives show a trend toward their 
respective mean value: 

Mt > p : Mi(p+1) = A(Ei -p)(m-p), Mi{p) = 1 - A(Ei -p){m-p); 
Ml <p: Mi{p-1) = X{p-E 2 ){m-p), Mu{‘2.,p) = l-\{p-Ei){m-p), 

where i, 2 . 

Therefore the model needs five phenomenological parameters, namely (xa 
is the susceptibility of citizens to become criminals; as is the susceptibility of 
criminals to reach back the state of normal citizens; /3 refers to the learning 
dynamics among criminals; 7 models the motivation/efHcacy of security forces 
to catch criminals; and A the learning (attraction) dynamics with respect to the 
mean value. 

Let us now briefly comment the result of simulations focusing on the role 
of the wealth distribution. The approach allows to investigate not only the role 
of macroscopic quantities such as the mean wealth, but also its distribution at 
equal mean value. More in detail, referring to the influence of the mean wealth 
the following result is obtained: 

• Decreasing mean wealth of the society Increasing the number of crim¬ 
inals and increasing their criminal ability; 

• Increasing mean wealth of the society Decreasing the number of crim¬ 
inals and decreasing their criminal ability; 

Let us now focus on the influence of the shape shape of wealth distribution 
the following result is obtained: 

• Equal distribution in a poor society Slow growth of the number of 
criminals; 

• Equal distribution in a wealthy society Fast decrease of the number 
of criminals; 

• Unequal distribution in a poor society Fast growth of the number of 
criminals; 

• Unequal distribution in a wealthy society Slow decrease of the number 
of criminals; 
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6.3 Possible developments 

As already mentioned, the model has been derived under quite limited assump¬ 
tions on the number of functional subsystems and in absence of a dynamics over 
the network. Both aspects of the dynamics are important as a higher number of 
levels of criminal behaviors and detective expertise are definitely closer to real¬ 
ity, while the effect of networks can have an important influence on the overall 
dynamics of the system. Indeed, social aspects of the interaction between wealth 
dynamics and criminal behaviors are object of intense studies by experts in the 
field of social science [72l UHl | 95 ] . 

We can observe that a number of recent papers, as examples P?lll 23 ] . study 
the concentration, and related properties, of hot-spots of criminals. A more de¬ 
tailed information can be achieved by subdividing the territory into a network of 
interconnected sub-territories and subsequently in studying the migration and 
concentration of criminals in certain zones. A challenging objective consists 
in deriving macroscopic scale equations from the underlying description at the 
microscopic scale. More in details, hyperbolic models with the feature of finite 
speed of propagation can be obtained [ 29 ] , as well as degenerate parabolic mod¬ 
els [ 5 D]. Previous applications to vehicular traffic m can be developed also in 
the case of social dynamics. 

However, according to the authors’ bias the study of the dynamics over 
networks is definitely more important also in view of its practical applications. 

As an example, if a certain territory, as an example a city, is subdivided in a 
number of subareas and if migration phenomena across them are modeled, the 
concentrations of criminality can be precisely identified and thus the planning 
of police forces distribution can be properly organized. 

7 Some Reasonings on the Validation of Models 

Focusing on a case study, it has been shown how the proposed approach can be 
applied to derive specific models. However, the general problem of model vali¬ 
dation has not yet found an exhaustive reply not only in the existing literature, 
but also in this present paper. Therefore, it is necessary tackling the last key 
question: 

How and how far mathematical models can be validated? 

Although not yet exhaustive, a partial answer will be given by going back to 
the key problems posed in [^. In general, validation requires that models can 
reproduce empirical data at the scale used to observe and represent the phe¬ 
nomena they describe. Moreover they should reproduce, at least at a qualitative 
level, emerging collective behaviors. 

The reasonings presented in the following do not claim to provide, an exhaus¬ 
tive reply on this challenging topic, but simply a focus on some key problems 
looking ahead to research perspectives. 


28 




i) Multiscale problems: Models might be designed at a scale different from 
that used for the observation and detection of empirical data. Generally, models 
refer to the microscopic scale, while only collective behaviors at the macroscopic 
scale are observed. 

ii) Dynamics of complex systems: Complex systems are characterized by 
collective behaviors which are subject to large deviations. However, qualita¬ 
tive shapes are in several cases preserved. Models are supposed to reproduce 
all observed emerging behaviors and, hopefully, also to predict, in some cases, 
behaviors that were previously unnoticed. 

iii) From interactions to collective dynamics: Collective dynamics are 
determined by interactions at the microscopic scale, which in turn are ruled by 
the strategy that interacting entities are able to express. Specific experiments 
should be designed to understand the dynamics of individual social interactions. 

iv) Evolution of rules of interactions with onset of rare events: A math¬ 
ematical approach to modeling should consider that individual behavioral rules 
and strategies are not, in most cases, constant in time due to the evolutionary 
characteristics of living complex systems. For instance, interactions among the 
individuals can change in time depending on several factors, including the global 
state of the system. Such variability can be a source of unpredictable events. 
Networks can enhance rapid spreads across nodes as a sort of domino effect 

Remark 7.1 Taleb’s definition U 29 \/ of the so-called “Black Swan” was already 
presented in the first section. It is worth stressing that this author is a rare exam¬ 
ple of researcher moving against the main stream of the traditional approaches, 
generally focused on well-predictable events. As such, it had an important impact 
on the guest for new research perspectives. In line with the approach followed 
in some speculations can be developed to provide an answer to the follow¬ 
ing key topic which deserves further speculations, that is the interplay among 
different dynamics an important ingredient for the onset of a non predictable 
event. 

All reasonings presented in this subsection do not directly contribute to the 
validation of models, but they enlighten the conceptual difficulties to be tackled. 
Indeed, the present state of the art does not yet provide well defined validation 
rules. Therefore, additional speculations can be presented by analyzing how far 
the theory proposed in this paper can provide a constructive answer to the key 
problems proposed in the book and reported in Section 2 . The reasonings 
proposed in the following correspond to the personal opinion of the authors of 
this present paper. 

1. Separation between homogeneous and heterogeneous behaviors: 

The separation between these two types of dynamics, which correspond 
to homophyly and heterophily, is modeled by a threshold, which sepa¬ 
rates consensus from dissent. According to the approach of this paper, 
the threshold depends on the distance between the interacting entities. In 
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some cases, for example opinion formation, small distances induce consen¬ 
sus, while the opposite trend is observed in the case of social competition. 
It is worth mentioning that such a threshold cannot be considered a con¬ 
stant quantity, as it depends not only on external actions, but also to the 
overall state of the system. 

2 . Dynamics toward consensus: Understanding the specific features that 
promote consensus has been a constant subject of research activity by 
mathematicians and physicists, examples are given by [Ml I 109 j . A recent 
survey [ 109 ] provides a deep analysis on this topic. The approach pre¬ 
sented in this paper links the dynamics toward consensus (cooperation) 
to the threshold already mentioned in the preceding item. Therefore, the 
trend increases (or decreases) according to the dynamics and external ac¬ 
tions, which modifies such a threshold. Papers [ 9 l [ 69 ] present specific 
applications related to this topic. 

3 . Dynamics over networks: The main observation concerning the model¬ 
ing of networks is related to the way by which the partition into functional 
subsystems has been organized. As already mentioned, this process is not 
developed according to universal rules, but it is related on the specific 
investigation which is pursued. As a consequence, the organization into 
functional subsystems and hence into a network is not fixed for all times 
as the dynamics of subsystems might induce changes in the network to be 
viewed as a dynamical system [102] . where new subsystems can appear, 
while others disappear. 

4 . Measure of centrality of networks: This topic remains open after the 
contribution of this present paper, which is limited to analyze how the 
network as a whole has an influence over the dynamics of each functional 
subsystem. A quantitative indication has been related to the influence 
of the mean value. However, this topic deserves additional studies to 
understand the role of higher order moments as well as the higher or 
lower influence of one node with respect to the others. This aspect can 
be accounted for by the interaction rate, while an interesting objective of 
investigation consists in understanding how far the action of the network 
can modify the dynamics within each node. 

5. Unpredictability and chaotic dynamics: This topic deserves atten¬ 
tion and should be developed according to the class of equations used in 
the modeling approach. Stability and chaotic dynamics are rapidly related 
to models stated in terms of ordinary differential equations. Focusing on 
the stochastic equations proposed in this paper, a natural investigation 
is the analysis, by appropriate entropy functions [ 18 ], of the asymptotic 
trend by understanding whether shows a trend toward a deterministic be¬ 
havior or toward an increasing variance. At present only computational 
results support this study [ 3 H [ 35 ] , while analytic studies definitely deserve 
future contributions. 
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8 Looking ahead to Research Perspectives 

This paper does not simply focus on a survey of the existing literature but it 
proposes a systems approach to social sciences. Such approach is developed in 
two steps. The first step consists in the derivation of a general mathematical 
structure suitable to capture the main features of behavioral social systems and 
the second step in deriving specific models by specializing the parameters of the 
structure according to each specific system under consideration. This approach 
includes the dynamics over networks and migration dynamics across the various 
nodes of the network. 

Looking ahead to research perspectives, firstly we point out the advantages of 
the aforementioned systems approach and then we critically analyze the shades 
to be hopefully enlightened by future research activity. 

According to the authors’ bias, the main advantage is the broad generality 
of the mathematical structure, which appears to be consistent with the main 
features of living systems. Moreover, it has been verified that the key features 
proposed in the book [ 53 ] have been effectively treated and a constructive answer 
has been given. Therefore, it can be stated that such a structure offers a tool 
useful for the applications as it has been shown in the case study reviewed in 
Subsection 5 . 2 . 

Possible improvements of the approach can be obtained by studying into 
more depth hve specific topics, selected according to the authors’ bias and are 
treated in the following. These topic share a common focus on the modeling 
interactions according to the idea that some preliminary reasonings on this topic 
can hopefully contribute to future research activity. 

8.1 Modeling interactions to-ward a socio-mathematical theory: A deep 
understanding of the dynamics of interactions can provide an important con¬ 
tribution to the study of collective dynamics. The mathematical structure pro¬ 
posed in this paper can transfer, as we have seen, the individual behavior to the 
collective dynamics. Theoretical contributions from social sciences can provide 
a theoretical approach to such a modeling. This challenging objective can estab¬ 
lish a constructive link between mathematics and social sciences by linking the 
mathematical structure to theories of social systems with some analogy to the 
interaction between physics and mathematics which, in the past century, gen¬ 
erated the so-called mathematical physics. The final objective consists in the 
development of a social-mathematical theory suitable to offer a vision broader 
than that offered by single models. 

8.2 Learning dynamics: This concept has been mentioned several times in 
the paper. Indeed, it plays a key role on the modeling of interactions, as active 
particles modify their rules of interactions according to what they learn from 
other particles and has been already object of attentions in various applications 
as we have seen in Subsection 5 . 2 . Accordingly, the theoretical approach to 
learning is not simply collecting and organizing a large set of information, as 
it needs modeling the actual dynamics by which individuals “learn” from the 
others. This approach has been recently developed in [ 56 ] based on the idea 
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that understanding the dynamics can contribute to modeling interactions. Par¬ 
ticularly important is the modeling of learning dynamics across networks [2], 
where individuals perceive the overall state of the functional subsystems and 
learn out of this knowledge with a direct influence over the social dynamics. 

8.3 Multiple dynamics: Some recent papers, see among others, [36l|69], have 
shown that the interactions of different dynamics can explain social phenomena 
that cannot be foreseen when each dynamics is taken into account separately 
from the other. Interesting examples are treated [351 US]. This literature mo¬ 
tivates further research in this field, where a challenging objective consists in 
understanding how the interaction rules of one dynamics can be extended to 
the case of multiple dynamics. Generally, a factorization of the transition prob¬ 
ability density is assumed to simplify mathematical models, but this heuristic 
assumption does not find any theoretical support. 

8.4 Additional developments on dynamics over networks: This paper 
has already proposed some perspective ideas on the modeling of networks by 
selecting various types on interactions either within the same node or across 
them. Two scales have been used, namely the microscopic one related to inter¬ 
actions involving active particles and micro-macro interactions between active 
particles and nodes as a whole. However, some limits still need to be prop¬ 
erly understood. Examples include the limitation to first order moments in the 
micro-macro interactions and the assumption, in the micro-interactions, that 
each particle interact with all particles in a node. The use of higher order mo¬ 
ments can contribute to model more precisely interactions. As an example the 
interaction rate could be modeled based on the conjecture that FSs featured by 
a high second order moment (somehow corresponding to energy) have a stronger 
influence on the FSs with low second order moment, namely the encounter rate 
grows with the difference between the relative moments. Moreover, interactions 
might be selective and limited to a low number of particles due to topological, 
rather than metric, interactions. This concept is proposed in m- 

8.5 Analytic problems: The application of models to the study of real prob¬ 
lems generates some challenging analytic problems, most of them have already 
been mentioned in the preceding sections focused on the initial value problem 
for specific models. We have already stressed that the proof of solutions to such 
a problem simply needs the application of standard techniques. Paper m has 
shown under which assumptions, reasonable from the view point of applications, 
existence and uniqueness can be proved by almost straightforward applications 
of fixed point theorems m- In fact, the structure of the equations guarantees 
local existence, under appropriate Lipschitz continuity assumptions of the in¬ 
tegral operators. Moreover, the solutions can be prolonged for arbitrary long 
times as the aforementioned structure assure preservation of the Li norm. 

On the other hand, the study of the asymptotic trend of the solutions is still 
waiting for a satisfactory answer. Numerical evidence suggests that for each ini¬ 
tial condition, solutions tend, asymptotically in time, to a configurations which 
appears to be unique and stable. However, this phenomenological observation 
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is not followed by a proof. A simple reasoning observes that this specific class of 
systems need a new definition of equilibrium. Indeed asymptotic configurations 
depend on the shape of the initial condition. Moreover, we stress that the study 
of systems subject to external actions is still looking for detailed analysis related 
to existence and control of solutions. 

Of course the five topics, presented above, do not claim to offer an exhaustive 
panorama, but simply promote new research initiatives as well as contributions 
of interested readers. It is not surprising that all of them focus on interactions. 
In fact, all different mathematical approaches to active particles, such as meth¬ 
ods of Brownian particles [118] or kinetic theory methods |115| , have to tackle 
the conceptual difficulties of modeling interactions among individuals. Which 
is the hrst research topic presented in this section. 
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